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■ A local current of particle density for scalar fields in curved background 

is constructed. The current depends on the choice of a two-point function. 
There is a choice that leads to local non-conservation of the current in a time- 
dependent gravitational background, which describes local particle production 
consistent with the usual global description based on the Bogoliubov transfer- 
ee \ mation. Another choice, which might be the most natural one, leads to the 
O" 1 local conservation of the current. 
U 
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One of the main problems regarding quantum field theory in curved space-time is 
how to introduce the concept of particles. The problem is related to the fact that 
particles are global objects in the conventional approach jl], |2|, while covariance with 
respect to general coordinate transformations requires local objects. This led some 
experts to believe that only local operators in field theory were really meaningful and, 
consequently, that the concept of particles did not have any fundamental meaning 
0. However, particles are what we observe in experiments. Moreover, from the 
experimental point of view, nothing is more local than the concept of particles. If 
we require that quantum field theory describes the observed objects, then it should 
describe particles as local objects. In this letter we show that the concept of particles 
can be introduced in a local and covariant manner. 

Our approach is based on a similarity between the number of particles and charge. 
For complex fields, the total number of particles is the sum of the number of particles 
and antiparticles, while the total charge is the difference of these two numbers. The 
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concept of charge can be described in a local and covariant manner because there 
exists a local vector current of charge density. We find that a similar vector current 
exists for the number of particles as well. Nevertheless, it appears that this local cur- 
rent is not unique, but depends on the choice of a two-point function W(x, x'). When 
a unique (or a preferred) vacuum |0) exists, then W(x,x') is equal to the Wightman 
function (0\<j)(x)<j)(x')\0) and the current is conserved. When such a vacuum does not 
exist, then there is a choice of W that leads to local non-conservation of the current in 
a time-dependent gravitational background, which describes local particle production 
consistent with the usual global description based on the Bogoliubov transformation. 
Another choice, which might be the most natural one, leads to the local conservation 
of the current in an arbitrary curved background, provided that other interactions 
are absent. 

Let us start with a hermitian scalar field in a curved background. The field 
satisfies the equation: 

(V^ + m 2 + 6R)0 = O, (1) 

where R is the curvature. Let £ be a spacelike Cauchy hypersurface. We define the 
scalar product 

' dSPfcdpfa, (2) 



s 

where 9 M is the usual antisymmetric derivative. If 0i an d 02 are solutions of flU), then 
(0) does not depend on S. It is convenient to choose coordinates (i, x) such that 
t — constant on S. In these coordinates, the canonical commutation relations can be 
written as [0(x), 4>(x')]x = [d (f)(x), & (j)(x')]-£ = and 

= dVz5 3 (x-x'). (3) 

The label £ denotes that x and x' lie on E. The tilde on = \g^\ 1 ^ 2 n tJ ' denotes 
that it is not a vector, while n M is a unit vector normal to E. 

Let us first construct the particle current assuming that a unique (or a preferred) 
vacuum exists. The field can be expanded as 

4>{x) = a kfk(x) + a\fl{x) , (4) 

k 

where (f k , f k ,) = -(/£, /£,) = S w , (/£, /*/) = {f k , /£) = 0. The operators 

flfc = (A, 0) , 4 = 4>) (5) 

satisfy the usual algebra of lowering and raising operators. This allows us to introduce 
the vacuum as a state with the property a k \0) = and the operator of the total 
number of particles as 

AT = £4a fc . (6) 

k 

We also introduce the function W(x,x'), here defined as 

W(x,x>) = Y.fk(x)r k (x'). (7) 

k 
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(Later we also study different definitions of W(x, x').) For future reference, we derive 
some properties of this function. For the fields described by we find 

W(x,x') = (O|0(x)0(x')|O) . (8) 

From fl7|) we find W*(x, x') = W(x',x). From the fact that f k and /£ satisfy Eq. (|I|) 
we find 

(V^d^ + m 2 + £R(x))W(x,x') = , (9) 

(V'^ + m 2 + ^ J R(x , ))W / (x,x / ) = 0. (10) 

From (|8|) and the canonical commutation relations we find that /j. and /£ are functions 
such that 

W(x,x')|s = WV,x)|e, 
9 ^(x,x')|e = d d' Q W(x',x)\v. (11) 

n°^[iy(x, x') - W(x', x)] s = «5 3 (x - x') . (12) 

So far nothing has been new. However, a new way of looking into the concept of 
particles emerges when (|5|) is put into (|6]) and (H) and (0) are used. This leads to the 
remarkable result that ([|) can be written in the covariant form as 

N = J^j^x) , (13) 

where the vector j^(x) is defined as 

j^x) = Jj^{W(x,x')d^ v cf>{x) ( f>{x , ) + h.c.}. (14) 

Obviously, the hermitian operator j^,{x) should be interpreted as the local current of 
particle density. If ([|) is used in (0), one can show that j M (a;)|0) = 0. 

Assuming that (||) is the usual plane-wave expansion in Minkowski space-time and 
that the (t, x) coordinates are the usual Lorentz coordinates, we find 



J M = f d 3 xj^(x) = ^2 — 4 a k • (15) 

J k W k 



The quantity k % /uj^ is the 3-velocity v\ so we find the relation J\n q ) = vn q \n q ), 
where \n q ) is the state with n q particles with the momentum q. The relations of this 
paragraph support the interpretation of j M as the particle current. 

Note that although j^{x) is a local operator, some non-local features of the particle 
concept still remain, because (|14"D involves an integration over S on which x lies. 
Since <f>(x') satisfies (^) and W(x,x') satisfies fllPp, this integral does not depend on 
E. However, it does depend on the choice of W(x,x'). Note also that the separation 
between x and x' in flT4]) is spacelike, which softens the non-local features because 



W(x, x') decreases rapidly with spacelike separation. As can be explicitly seen with 
the usual plane-wave modes in Minkowski space-time, W(x, x') is negligible when the 
spacelike separation is much larger than the Compton wavelength m -1 . 
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Using ([]]) and (jfty, we find that the current ( |T4| ) possesses another remarkable 
property: 

V%(x) = . (16) 

This covariant conservation law means that background gravitational field does not 
produce particles, provided that a unique (or a preferred) vacuum exists. 

To explore the analogy with the current of charge, we generalize the analysis to 
complex fields, with the expansion 

<t>i x ) = ^ a kfk{x) + b\fl(x) , 

k 

4>Kx) = Y.a\r k {x) + b k f k {x) . (17) 

k 

We introduce two global quantities: 

N (±) = £ 4a* ± (18) 
k 

Here is the total number of particles, while is the total charge. In a similar 
way we find the covariant expression 

= I dW 3 ^\x) , (19) 



where 

jf)(x) = jj^\{W(x,x')d,d' v 

[<P\x)<P(x') ± <p(x)<p\x')} + h.c.}. (20) 

The operator jj^'(x) is the generalization of (|I~4"D with similar properties, including 
the non-local features. On the other hand, the apparent non-local features of jk~'(x) 
really do not exist, because, by using the canonical commutation relations, fll!]) and 
(|j~2|), j^'(x) can be written as 

jjr\ x ) = i(f>\x)d fl (j)(x) 

+ dY! v W[x, x') dj? v W(x', x) , (21) 

so all non-local features are contained in the second term that does not depend on 
<p. Using (0), (0) and the orthonormality relations among the modes f k , one can 

show that the quantity — (O|i0T(x)c^0(a;)|O) is equal to the second term, which reveals 
that the second term represents the subtraction of the infinite vacuum value of the 
operator represented by the first term. In other words, is the usual normal 
ordered operator of the charge current. 

In a way similar to the case of hermitian field, we find the properties jy~'\0) = 
and V^j^ = 0, while the generalization of (|l^) is 

= I d * x ^ (*) = E — (4«k ± blb k ) . (22) 



k 
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Let us now study the case in which a non-gravitational interaction is also present. 
In this case, the equation of motion is 



(V M <9 M + m 2 + ^R)4> = J , (23) 

where J(x) is a local operator containing (f> and/or other dynamical quantum fields. 
Since it describes the interaction, it does not contain terms linear in quantum fields. 
(For example, J(x) may be the self-interaction operator — A0 3 (s)). We propose that 
even in this general case the currents of particles and of charge are given by the 



expressions (|14j) and (P0|) , where W(x,x') is the same function as before, satisfying 
the "free" equations (||) and (|10|) . As we show below, such an ansatz leads to particle 
production consistent with the conventional approach to particle production caused 
by a non-gravitational interaction. 

Note also that our ansatz for W makes (21) correct even in the case J ^ 0, because 



the interaction does not modify the canonical commutation relations. This implies 
that the charge is always conserved, provided that the Lagrangian possesses a global 
U(l) symmetry. On the other hand, using (0) and fl2"3]), we find 



= j^<m' vl -{W{x,x')& v J(x)4>(x') + h.c.} , (24) 

and similarly for Note that only x (not x') appears as the argument of J on 

the right-hand side of (|24f) , implying that J plays a strictly local role in particle 
production. 

Let us now show that our covariant description (|4|) of particle production is 
consistent with the conventional approach to particle production caused by a non- 
gravitational interaction. Let £(i) denote some foliation of space-time into Cauchy 
spacelike hypersurfaces. The total mean number of particles at the time t in a state 
|V>) is 

tf(t) = (1>\N(t)\4>) , (25) 

where 

N(t)= [ dE*%(t,x). (26) 

Equation ( |25| ) is written in the Heisenberg picture. However, matrix elements do 
not depend on picture. We introduce the interaction picture, where the interaction 
Hamiltonian is the part of the Hamiltonian that generates the right-hand side of (|2"3|). 
The state \if)) becomes a time-dependent state \ip(t)), the time evolution of which is 
determined by the interaction Hamiltonian. In this picture the field <fi satisfies the 
free equation ([]]), so the expansion (f|) can be used. Since we have proposed that W 
satisfies the free equations @ and fllCf ), this implies that N(t) becomes J2k a l a k m 
the interaction picture. Therefore, (p5|) can be written as 



W) = ZW)\a]flk\m), (27) 

k 

which is the usual formula that describes particle production caused by a quantum 
non-gravitational interaction. 
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Let us now show that ( |HD can also describe particle production by a classical 
time-dependent curved background, provided that we take a different choice for the 
function W(x,x'). When the particle production is described by a Bogoliubov trans- 
formation, then the preferred modes fk in (f|) do not exist. Instead, one introduces 
a new set of functions ui(x) for each time t, such that ui(x) are positive- frequency 
modes at that time. This means that the modes ui possess an extra time dependence, 
i.e. they become functions of the form ui(x;t). These functions do not satisfy ([]]). 
However, the functions ui(x;r) satisfy ([1]), provided that r is kept fixed when the 
derivative <9 M acts on u\. To describe the local particle production, we take 

W{x, x') = Mx; t)u*{x'; t') , (28) 
i 

instead of (0). Since ui(x;t) do not satisfy ([]]), the function ([^) does not satisfy (|9|). 
Instead, we have 

(V M <^ + m 2 + £R(x))W(x, x') = -K(x, x') ^ . (29) 
Using ([!]) and (|29D in fll4]) , we find a relation similar to (|24l): 

V%(x) = J d^ ,v ^{K{x,x')& v (j>{x)(j>{x') +h.c.} . (30) 

This local description of the particle production is consistent with the usual global 
description based on the Bogoliubov transformation. This is because fl2"E| ) and ( |HD 
with (|28|) and (i) lead to 



N(t) = ^A}(t)Mt), (31) 
i 

where 

Mt)=J2^ k (t)a k -P: k (t)a[ 1 (32) 

k 

<*ik(t) = (fk,ui) , Pik(t) = -(fk,ui) . (33) 

The time dependence of the Bogoliubov coefficients an~(t) and $&(£) is related to the 
extra time dependence of the modes Ui(x;t). If we assume that the change of the 
average number of particles is slow, i.e. that d t Ai(t) ~ 0, d t Ui(r, x; t)\ T=t ~ 0, then 
the Bogoliubov coefficients ( |3"BD are equal to the usual Bogoliubov coefficients. This 
approximation is nothing else but the adiabatic approximation, which is a usual part 
of the convential description of particle production Q. 

In general, there is no universal natural choice for the modes ui(x; t). In particular, 
in a given space-time, the choice of the natural modes ui(x;t) may depend on the 
observer. If different observers (that use different coordinates) use different modes for 
the choice of (p8|) , then the coordinate transformation alone does not describe how 
the particle current is seen by different observers. In this sense, the particle current 
is not really covariant. There are also other problems related to the case in which 
different observers use different modes [Q. 

Since the covariance was our original aim, it is desirable to find a universal natural 
choice of W(x, x'), such that, in Minkowski space-time, it reduces to the usual plane- 
wave expansion (|7|). Such a choice exists. This is the function G + (x, x') satisfying 
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(0) and (TO), calculated in a well-known way from the Feynman propagator G F (x, x') 
0. |j. The Feynman propagator G F (x,x'), calculated using the Schwinger-DeWitt 
method, is unique, provided that a geodesic connecting x and a;' is chosen. When x 
and x' are sufficiently close to one another, then there is only one such geodesic. In 
this case, the adiabatic expansion |2], |6| of the Feynman propagator can be used. For 
practical calculations, the adiabatic expansion may be sufficient because G + (x,x') 
decreases rapidly with spacelike separation, so the contributions from large spacelike 
separations may be negligible. 

To define the exact unique particle current, we need a natural generalization of 
G F (x,x') to the case with more than one geodesic connecting x and x'. The most 
natural choice is the two-point function G F (x,x') defined as the average over all 
geodesies connecting x and x' . Assuming that there are N such geodesies, this average 
is 

N 

G F {x, x') = N- 1 ]T G F {x, x'- a a ) , (34) 

a=l 

where G F (x, x'\ a a ) is the Feynman propagator G F (x,x') calculated with respect to 
the geodesic a a . Eq. (^4j) can be generalized even to the case with a continuous 
set of geodesies connecting x and x' . This involves a technical subtlety related to 
the definition of measure on the set of all geodesies, which we shall explain in detail 
elsewhere. 

Since the function G + (x, x') calculated from G F (x, x') always satisfies @ and (|T0|), 
it follows that the corresponding particle current is conserved and does not depend 
on the choice of E, provided that the field satisfies ([!]). Note that this definition of 
particles does not always correspond to the quantities detected by "particle detectors" 
of the Unruh-DeWitt type 0, § . Instead, this number of particles is determined by a 
well-defined hermitian operator that does not require a model of a particle detector, 
just as is the case for all other observables in quantum mechanics. Moreover, since ( |T4"D 
does not require a choice of representation of field algebra, the definition of particles 
based on Gf does not require the choice of representation either. This allows us 
to treat the particles in the framework of algebraic quantum field theory in curved 
background |J. Furthermore, as noted by Unruh |7[], only one definition of particles 
can correspond to the real world, in the sense that their stress-energy contributes 
to the gravitational field. Since the definition of particles based on G F is universal, 
unique and really covariant, it might be that these are the particles that correspond 
to the real world. Besides, as we shall discuss elsewhere, it seems that these particles 
might correspond to the objects detected by real detectors (such as a Wilson chamber 
or a Geiger-Miiller counter) in real experiments. 

In this letter we have constructed the operator describing the local current of 
particle density, which allows us to treat the concept of particles in quantum field 
theory in a local and generally covariant manner. This operator is not unique, but 
depends on the choice of the two-point function W(x,x'). Different choices corre- 
spond to different definitions of particles. In particular, various choices based on ( p8|) 
correspond to particle production by the gravitational field. This local description 
of particle production is consistent with the conventional global description based on 
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the Bogoliubov transformation. Similarly, various choices based on (|7|) give a local 
description of particle content in various inequivalent representations of field algebra. 

We have seen that a particularly interesting choice of W(x, x') is that based on 
the Feynman propagator Gf(x, x'), because this choice might correspond to the most 
natural universal definition of particles for quantum field theory in curved background. 
It is tempting to interpret these particles as real particles. If this interpretation is 
correct, then classical gravitational backgrounds do not produce real particles. There 
are also other indications that classical gravitational backgrounds might not produce 
particles flO], [□], [T^|. Even if the particles based on Gf do not correspond to real 
physical particles in general, it is interesting to ask about the physical meaning of 
this time-independent hermitian observable that corresponds to physical particles at 
least in Minkowski space-time. 

Note finally that our definition of the particle current can be generalized to the 
case of scalar and spinor fields in classical electromagnetic backgrounds. Again, a 
two-point function can be chosen such that the particle production is described in a 
way consistent with the Bogoliubov-transformation method, but the requirement of 
gauge invariance leads to a conserved current, indicating that classical electromagnetic 
backgrounds also might not produce particles. Only the quantized electromagnetic 
field (having a role similar to J(x) in ([23])) may cause particle production in this 
formalism, which is in agreement with some other results [Q. This will be discussed 
in more detail elsewhere. 

This work was supported by the Ministry of Science and Technology of the Re- 
public of Croatia under Contract No. 00980102. 
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